Quantum radiation in a plane cavity with moving mirrors 
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We consider the electromagnètic vacuum field inside a perfect plane cavity with moving mirrors, 
in the nonrelativistic approximation. We show that low frequency photons are generated in pairs 
that satisfy simple properties associated to the plane geometry. We calculate the photon generation 
rates for each polarization as functions of the mechanical frequency by two independent methods: 
on one hand from the analysis of the boundary conditions for moving mirrors and with the aid of 
Green functions; and on the other hand by an effective Hamiltonian approach. The angular and 
frequency spectra are discrete, and emission rates for each allowed angular direction are obtained. 
We discuss the dependence of the generation rates on the cavity length and show that the effect is 
enhanced for short cavity lengths. We also compute the dissipative force on the moving mirrors and 
show that it is related to the total radiated energy as predicted by energy conservation. 

os: 

^ . 42.50.Lc, 42.50.Dv, 03.65.-w 

< 
OC 

■ I. INTRODUCTION 

r-H ' 

' In the presence of moving boundaries, the vacuum state of the electromagnètic field may not be stable, then 
resulting in the generation of photons. This purely quantum effect, which has been known either as dynamical 
Casimir effect [|p or as motion or mirror || induced radiation is, like the usual Casimir effect for standing mirrors, 
a striking illustration of the physical reality of the quantum vacuum field. Moreover, it may also be understood 
as a mechanical effect of the vacuum field. In fact, energy conservation entails that the radiation effect must be 
accompanied by a radiation reaction force that works against the motion of the mirror j6|, and which is connected 
. to the fluctuations of the usual (static) Casimir force by the fluctuation-dissipation thcorem |Q- ||. 

Several theoretical models have been analyzed. In the one-dimensional approximation (1D), only one direction of 
propagation is taken into account JÏOf . The quantum radiation generated inside a 1D cavity with moving mirrors was 
calculated in Refs. [ïïf| and [ï^| in the particular case where the mechanical frequency satisfies a resonant condition for 
generation of photons in the lowest order cavity field modes, whereas Ref. Q considered a 1D cavity with partially- 
transmitting mirrors and with no particular assumption about resonance, thereby allowing for a full analysis of the 
^ | spectrum of the radiation in a more general case. 

A few three-dimensional (3D) models have been recently analyzed in the literature, including moving dielectric half- 
spaces |Q, and rotating |Ï5| or collapsing dielectric spheres |Ï6j, the latter as a model for sonoluminescence. On 
the other hand, 3D results for the related problem of photon generation in a medium with time-dependent material 
coefficients e and \i have been known for nearly ten years |Ï7j . Perhaps the simplest 3D illustration of motion induced 
radiation is to consider a single perfectly-reflecting plane mirror moving in free space. In the perturbative regime, 
which is associated to the nonrelativistic limit, it is possible to derive simple results for the spectra of radiation 
JÏ8| , which display interesting polarization-dependent features connected to the angular distribution of the emitted 
photons. In this paper, we extend the method developed in Ref. |ïg[| to analyze the radiation emitted when two 
parallel plane perfectly reflecting mirrors, initially a distance L apart, oscillate along the direction perpendicular to 
their surfaces, and according to a pre-defined law imposed by some external apparatus. Such geometry constitutes 
the simplest example, from a theoretical point-of-view, of a 3D cavity of length L. As compared to the previous 
single-mirror case, we show that the orders of magnitude for the radiation rates generated in the plane cavity may 
be several orders of magnitude larger, provided that L is small enough. 

The paper is organized as follows. In Sec. II, we calculate the photon numbers generated inside the cavity starting 
from the boundary conditions of a moving perfectly reflecting mirror. The method is based on the nonrelativistic and 
long-wavelength approximations, which are closely connected in the context considered here Q. In Sec. III we present 
an alternative derivation of the results already found in Sec. II, now employing usual time-dependent perturbation 
thcory for an effective Hamiltonian which incorporates the motion effect in terms of a coupling via radiation pressure. 
This heuristic approach is considerably simpler than the previous one, since it circumvents the analysis of the moving 
boundaries. Furthermore, it explicitly unveils the two-photon nature of the photon emission process, and allows for 
the computation of the dissipative component of the radiation pressure force on the moving mirrors. In Sec. IV, we 
consider an specific example of motion in order to isolate the effect of a single mechanical frequency lüq. We show that 
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the photon numbers obtained by two independent methods in Secs. II and III grow linearly in time, allowing us to 
define photon production rates, whose behavior as functions of the dimensionless parameter lüqL/ttc is examined in 
detail. Sec. V contains the concluding remarks. 



II. BOUNDARY CONDITIONS AND INTRACAVITY QUANTUM RADIATION 

For the sake of clarity we first assume that one of the mirrors is at rest. The results in the more general case where 
both mirrors are set to move is a simple generalization to be presented later. The moving mirror oscillates along the 
direction perpendicular to its surface (x direction), around the position x = 0, its instantaneous position being given 
by the equation x = Sq(t). 

We decompose the electromagnètic fields into their components corresponding to the elèctric field parallel (TM) or 
perpendicular (TE) to the plane of incidence. For each polarization it is possible to define a vector potential through 
the equations: 

E (TE) = _Q tA (T E) . B (TE) = y x A (TE) ^ 

and 

E (TM) = y x ^(TM). B (TM) = ^(TM)^ ^ 

The units are MKS with c = 1 and eq = 1. The potentials satisfy the Gauge equations 

V • A (TE) = V • A iTM) = 0. (3) 

As shown in Appendix A, the boundary conditions for a perfectly reflecting moving mirror are very simple when 
written in terms of A (TE) and *4 (TM) , due essentially to the fact that they are both orthogonal to the direction of 
motion. We find 

A^(x = Sq(t),r ll· t)=0 (4) 

and 

(d x + Sq(t)dt)A (TM \x = Sq(t),r h t) = 0, (5) 

where r || = yy+zz. Furthermore, the fields satisfy the usual homogeneous Dirichlet and Neumann boundary conditions 
on the second mirror, which is at rest at x = L : 

A< TE > (x = L, r,| , í) = 0; 8 X A TM) (x = L,r^,t) = 0. (6) 

We want to solve the boundary value problem as defined by Eqs. (§)-(||) for the fields in the region between the 
mirrors. The results for the fields outside the plane cavity are essentially the same as those for a single moving 
mirror in vacuum, and hence may be found in Refs. |q,p^|. The essential 'ansatz' which allows us to employ the long 
wavelength approximation to solve the boundary value problem defined by Eqs. (||)~(^|) is to assume that a given 
mechanical frequency wo induces the generation of photons only in the spectral range ui < lü . This property is satisfied 
by the nonrelativistic models considered previously (see Refs. ^ and |Ï8|]). Moreover, it agrees with the intuitive 
notion that the radiation effect is a nonadiabatic process, so that high frequency field modes cannot be excited since 
the corresponding time scales are shorter than mechanical time scales (quasi-static limit). More importantly, we show 
later in this section that this property is fully satisfied for the model considered here. As for the connection with the 
long wavelength approximation, we note that the amplitude Sqo of a sinusoidal nonrelativistic motion must satisfy 
LüoSqo <C 1. When combined with our ansatz, this condition leads to <5go A, where A is the wavelength of the emitted 
radiation. Actually, we may be slightly more general and consider any nonrelativistic oscillatory motion around x = 
such that its Fourier components satisfy the above requirements (more specifically, we shall consider a weakly-damped 
sinusoidal nonrelativistic motion in Sec. IV). 

Accordingly, we look for perturbative solutions in the form: 

A< TE > = A^+5A^\ (7) 

and 
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A (TM) =Ag^ +SA <Tm . (8) 

Ag^a ' an d «^-s™' are the fields satisfying the Dirichlet and Neumann boundary conditions for standing mirrors, whereas 
5A {TE) and à4. ( ™ > represent the first-order modifications induced by the motion. As we show below, they are smaller 
than the fields for the static configuration by a factor of the order of Sq/X. We expand the fields in Eqs. (|4|) and (|5|) 
in Taylor series around x = 0. Since the j-th spatial derivative of a monochromatic travelling wave satisfies 

\díA\ < (Çy'|A|, (9) 

we find from Eq. (ji|) that the TE-polarized field SA (TE) is given up to first order in Sq/X by 

<5A< TE >(* = 0,r||,í) = -Sq(t)d x A^(x = 0,r h t). (10) 

Note that we have neglected the term Sq(t)d x SA ITE) (x — 0,rii,t) because, as shown by the above result, <5A (TE) is 
already of first order in Sq/X. Following the same method we find the following result for TM polarization: 

d x 5A™(x = 0,r h t) = -(Sq(t)d 2 x + Sq(t)dt)Ai™ } (x = O.r,,,*), (H) 

where now we have also neglected terms of the order of SqSq/X. According to our 'ansatz', when considcring the 
generation of photons out of the vacuum field induced by a mechanical frequency lüq, the relevant wavelengths 
are larger than 2ir/uio, and thus the neglected terms are all of the order of (Sq) 2 . We have then transformed the 
homogeneous boundary conditions for the total fields at the time dependent position x = Sq(t) given by Eqs. (|J) and 
(||) into inhomogeneous boundary conditions for <5A (TE) and à4. <TM) at the position x — 0, given byEqs. © and @, 
which may be solved by Standard Green function techniques. 

We introduce periòdic boundary conditions on the transverse plane yz over a surface of area S. In the static case, 
the normal mode decomposition of the fields in the interval between the mirrors, < x < L, is then written as follows: 




Ajr (.'./. = ,y y J-^sin^e"^ - '„,/: <„ II,. (I2! 



and for the TM polarization, 



-4£"M) ='EE^ (U fa * K5£ <M^'·e-"'·<^ + H.C (13) 

where 

kjf = 2n{n v y + n z z)/VS (14) 

represents the component of the wavevector parallel to the mirrors — the shorthand n — (n y , n z ) representing a pair 
of integer numbers. Note that the two potentials describing orthogonal polarizations are written in terms of the same 
unit vector 

k n 

èn=*X-I (15) 



Throughout the paper, the sum over n — as in Eqs. (12) and (|13j) — runs from n y — — oo and n z = — oo to n y = oo 
and n z — oo. A given mode with indexes (n, i) corresponds to a standing wave along the x-direction with wavevector 
k x — l-K jL travelling along a direction parallel to the mirrors with wavevector kí?. Its frequency is given by 



tir\ 2 (2?r) 2 
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The bosonic field operators in Eqs.([ï^) and (^3|) satisfy the usual commutation relations 



(%) 2 + [n z f . (16) 



and 



<íX'.í'] = °. ( 17 ) 
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l í,V ( a nV') f ] = ^n,n"5/,/"5j,i' 



(18) 



where j = TE, TM represents the polarization. 

It is convenient to work with a mixed Fourier representation defined as 



| / dtj d 2 r||e~ ïk r r »e í " t A (TE, (a;,r||,í) 



(19) 



with an analogous expression for TM polarization. The Fourier transformed fields representing the motion induced 
perturbation satisfy the 1D Klein-Gordon equation 



( 



(d 2 x +co 2 ~(k^)SA^[x,u;]=0, 
dl +lü 2 - (k\!) 2 ) SA ( ™ } [x,uj} = 0, 



(20) 
(21) 



and the boundary conditions at x — and x = L given by Eqs. (|^), ( |Ï0| ) and (pr]). The resulting boundary value 
problem for TE polarization is solved with the aid of the appropriatc Dirichlet Green function: 



qD . y sin(^)sin(^) 



(22) 



where the plus (minus) sign in Eq. (p2) provides the retarded (advanced) Green function. The fields with TM 
polarization are obtained from the Ncumann Green function: 



(l + ^o)((w±ie) 2 -^ ) 



(23) 



We assume that the mirror moves during a finite time interval, then returning to its initial position at x = 0. As 



a consequence, we may define input and output fields, Ai, 



and A„ 



(TE) 



corresponding to the limit vàlues of very 



small and very large times (and likewise in the case of TM polarization), which satisfy the boundary conditions for 
a mirror at rest at x = 0. They are connected by a suitable combination of retarded (superscript R) and advanced 
(superscript A) Green functions: 



a. < TE >r 



> [x, w] + SA^ [x' = 0, u] [d x .GÏ*{x, x' = 0) 
-d x ,G°> A (x,x' = 0) . 

The TM output field Amtíí™' is related to the TM input field An^™' by a similar expression: 



(24) 



0,Lü] 



G^ R (x,x'=Q) 



A aa rfr ) [x,«;] = A^T'M ~ 8*6 
-G^ A (x,x' = 0) . 

From Eqs. (|||) and (|||) we find 

d x ,G^ R {x,x' = 0) - d x ,G%< A (x,x> = 0) = Ei=i 4sin(^) 

and 



(25) 



(26) 



G^(x,x' = Q)-G^ A (x,x' = 0) 



_ 27TÏ V^°° 

L L·t=Q (l+S m )u e n 



(27) 



In general, there are no monochromatic solutions for the problem of moving boundaries, and hence it is not possible to 



write a normal mode decomposition for the field in this case. However, since Ai n ™ and A Q 



satisfy the boundary 



conditions for two standing mirrors at x = and x = L, we may write their normal mode decompositions as in 
Eqs. (|Ï2j)-(|Ï3|), in terms of input and output bosonic operators ain™ and a ou t^ TE) (at this point our method is quite 
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similar to the approach developed in Ref. |17]] for the problem of time-dependent material coeficients). We then take 
the Fourier transform of Eq. (|l(]) and replace the result, jointly with Eq. ((2^), into Eq. (Q) in order to find the linear 
transformation between the input and output TE bosonic operators: 



Oout 



(TE) _ (TE) . ]_ \ 

ni — " m ni "T" t / , 



í-k l'ir 
L L 



n n 



■: r ' t 'l (te) 



(28) 



where Sq[to\ is the Fourier transform of Sq(t). The relation between TM operators is derived from Eqs. ([ïï|), ( p5|) and 
(p7|) in a similar way: 



(TM) (TM) 

a out„e — Qiiw 



l£[(l + í, )(l + ^o) 



£'=0 



(29) 



7 



r r ^ ^'l (TM) . v "|l ' ' n 71 s: \ l , í'm (tm)n + 



From Eqs. ( pg ) and ( j29| ) we may readily derive the number of photons generated inside the cavity as a quantum 
cffcct of the mirror's motion. As discussed below, the effect is associated to the creation operators appearing in the 



r.-h.-s. of Eqs. (28) and (|29|). We assume that the field is initially in the vacuum state. The motion of the mirror 
then excites a given number of photons N 3 n t with indexes n, í and polarization j. f is given by the corresponding 
output number operator averaged over the input vacuum state: 

Nl, = (0, in | (a ou 4)+a ou 4 | 0, in), (30) 

Replacing Eqs. ( |28| ) and (|2|) into ([30j) provides the photon numbers for each polarization: 



'=1 



vZ(t) (t) ròr'^ + ·í" , · (31) 



and 



, ^ («) 2 +c^'Y 



Since the freqüències uj^ are positive, we infer from Eqs. ( |3Ïj ) and ( |32| ) that a given mechanical frequeney ojq generates 
photons with freqüències o>^ < u>o, thereby justifying the ansatz employed in this section. 

From the above results we may directly calculate the photon production rates and then estimate the order of 
magnitude of the quantum radiation effect. Before addressing this question, however, we present a second derivation 
of Eqs. (En!) and (|32|), which is based on usual time-dependent Hamiltonian perturbation theory. Note that the 



invariance of the r.-h.-s. of Eqs. ( |3l| ) and (32) with respect to the permutation of i and £' suggests that the photons 
are emitted in pairs. That this is indeed the case is more clearly shown by this alternative approach, to be presented 
in the next section. 



III. CONNECTION WITH RADIATION PRESSURE 



Rather than considering the boundary conditions of a moving mirror, we follow in this section the heurístic approach, 
first presented in Ref. B, in which the effect of the mirror's motion is modelled by taking the perturbation Hamiltonian 

SH = -Sq{t)F, (33) 

where F is the field quantum operator representing the force on the moving mirror. Accordingly, SH corresponds to 
the energy supplied to the field by means of the vacuum radiation pressure effect. The total Hamiltonian of the field 
is 
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H = H {a) + 5H, (34) 

where the unperturbed Hamiltonian is written in terms of the bosonic field operators for a standing mirror (see 
Eqs. @ and (P)) as 

^ (0) =E E + (35) 

n,£ j=TE,TM 

As discussed elsewhere |Ï(J , a Hamiltonian approach is not rigorously consistent with the model of perfect reflectiveness 
considered here. However, this model may be considered as an approximation for dielectric mirrors with large refraction 
index n — for which a rigorous Hamiltonian model is available although such correspondence is not yet settled 
(according to Ref. ]Ï3|] , some unexpected results show up when taking the limit of large n). In any case, the formalism 
presented in this section is justified by comparing the results it provides with those obtained in Sec. II. 

The force operator is the integral over the surface of the mirror (at its rest position at x = 0) of the xx component 
of the Maxwell stress tensor: 

F = \j d 2 r||(S x (0+) 2 - J B||(0+) 2 ), (36) 

where the limit x — > is taken from positive vàlues of x as indicated above (as in the previous section, we do not 
analyze the effect of the field outside the plane cavity). Since F is a quadratic operator on the field, the perturbation 
Hamiltonian SH excites pairs of photons like in the problem of parametric amplification by a nonlincar medium. 
Thus, we consider a perturbed field state of the form 

l*>= E Hnej,n'ef}(t)\{nij,n'e'f}) + b(t)\0), (37) 

where we sum over all two-photon states \{n£j,n'£'j'}) (the symbols j and f representing the polarizations of the 
photons in a given pair {n£j, n'£'j'}). Note that each pair {n£j, n'£'j'} is included only once in Eq. (|37|), regardless of 
the ordering of the indices. 

We assume that at t — > — oo the field is in the vacuum state, so that the two-photon amplitudes are initially 
zero: C{ n £j in >eijiy(— oo) = 0, b{— oo) = 1. We compute the build-up of the two-photon amplitude ci n tj ín >i'j'y(t) from 
Standard first-order perturbation theory: 

/* 
(n£j, n'£'f\SH(t')\0) exp - E^ c )t'\ dt', (38) 

with 

E$ n «,-EW=h(u>' n + u£,) (39) 

representing the difference between the (unperturbed) energies of the final and initial states. As discussed in the 
previous sections, it is mea ning less to discuss two-photon amplitudes as long as the mirror is moving. Accordingly, 
we must take t — > oo in Eq. (|3§|) in order to have a consistent picture of the quantum radiation effect. Then, replacing 
Eqs. (||) and (||) into Eq. (|§) yields 

c^n't'M 00 ) = ^jye'fmSqiüji + ü#]. (40) 

In order to compute the matrix element appearing in the r.-h.-s. of Eq. (pïo|), we write the elèctric and magnètic 
fields in Eq. ( |36] ) in terms of the potentials A (TE) and ^4. (TM) . It is convenient to use the Fourier series representation 
defined by 

A< TE »(x,r||,í) = ^A^(M)exp(ikf • r,,), (41) 

n 

and by an equivalent expression for the TM potential *45 TM) . Then, the force operator is written as 

F = f E [( fc ") 2 ^ TM) (o + , t) ■ ALT (o + , é) 

n 

-^ TM) (0 + ,t) ■d t A ( ™ ) (0 + 1 t) -d x A£ E) (0 + ,t) ■d x A^ l \0 + ,t)\. (42) 
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From Eq. (f42|), we obtain 

(nÍTE,n'£'TM\F\0) =0. (43) 

Therefore, the photons belonging to a given emitted pair have the same polarization. This is a general property of 
the plane symmetry of the problem, rather than a consequence of the specific model considered in this paper. Note 
however that it has been recently shown that TE-TM pairs may be radiated in the case of lateral motion of the 
mirror [fl4j| . 

Using the normal mode decomposition of the field operators as given by Eqs. (|ï^) and (|Ï3|), we may calculate the 
TE-TE and TM-TM matrix elements. We first find 

(n£TE,nr T E|^A^ E, (0+ ! í) .d x A™(0+,t')\0) = 

■fc ÍTX i' IX 



L L 



and 



V ^n^n 

(n€TM,nYTM|^ M) (0 + ,í)·^™ ) (0 + ^')|0) = ^(l + fco)(l + fco) 

x^X' (à n>N 6 n ,^ N e i ^+<- 1 : O + 5 n ^ N 8 nl , N e i ^ t+ < t '^ , (45) 

where use was made of the property iú l n — wl n . Combining Eqs. (Q) and ( pÏ5"|) with Eq. ( [42| ) leads to 

(nlT E ,nYTE|F|0) = -- ^_£_ J n ,_ n ,, (46) 



and 



fi (e) 2 + uiui 

(ní tm, tm | F 1 0) = " 5 n _„/ . (47) 

1(1+^(1 + ^0)^ 

From Eqs. ( j46|) a nd (|Í7|) , we may immediately calculate the amplitudes of creation of pairs of photons by combining 
them with Eq. ([40|). Here we write the results obtained in Appendix B for the more general case where both mirrors 
are moving, so that the first mirror is at x = <5çi(t) and the second mirror at x = L + «SfeCO- The resulting creation 
probabilities are: 

|2 1 (ÍTiV ÍÍ'tt\ 2 1 
|c Mte ,„^ te} | = J ? { T {— (48) 



n n 



and 



1 ((iqr + <J n <4) 



2 



\c {ntru , n >t>TH}\ L2 (1 + ^ q)(1 + 5lto]LJ t nUJ e x ( 49 ) 

Note that the photons in a given pair have opposite vàlues of , which is again a consequence of the plane symme- 
try p8| . As shown in Appendix C, Eqs. (H)~(^9|) must be slightly modified when considering the particular value 
n = n' = (which corresponds to the 1D limit of our 3D formalism, since such modes propagate along the x direction 
and do not contain any dependence on the transverse coordinates y and z). 

According to Eqs. ( fí8"| ) and (f49f), the joint motion of the two mirrors selects the longitudinal modes according to 
the parity of the indices í. When òq\ = —&[2, which corresponds to the 'elongation mode' of the cavity, the two 
photons in a pair correspond to l vàlues of the same parity, the opposite taking place when the motion is such that 
the cavity length is kept constant {bq\ = 5q2)- This property is a straightforward generalization of the situation found 
in one-dimensional cavities Q . It shows that the radiation effect should not be interpreted simply as a consequence 
of changing the optical cavity length, since it also takes place when there is no relative motion of the mirrors. 
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We may compute the average number of photons in a given cavity mode from 

<, = <*ltò)Vj*>. (50) 



Replacing Eq. (|37|) into Eq. (|50j) yields 

K,l = ^\c{nlj,-nt> 3 }\ 2 - (51) 



Eqs. ( f48| ) and ([49[), in thc particular case of Sq2 — 0, jointly with Eq. ( J5l| ) provide results for the photon numbers in 
full agreement with Eqs. (|3Ï|) and ( |3^ ) of the previous section. As for the particular case with n = 0, Eq. ( |5l"| ) also 
needs some slight modification in order to include the contribution of the degenerate two-photon states \0£,0£). As 
shown in Appendix C, there is agreement with the results found in section 2 in this case as well. We then conclude 
that the heuristic approach developed in this section yields the same final expressions for the number of photons 
produced in a given cavity mode. Moreover, it explicitly shows that the photons are generated in pairs, the photons 
in a pair having the same polarization and opposite vàlues of ky . 

With the aid of the linear response formalism |ï^] , the perturbation Hamiltonian as given by Eq. (|3^) may be also 
applied to compute the dissipative part of the radiation pressure force SF exerted on the moving mirrors ||- 
Such dissipative force is the mechanical effect of the quantum radiation process, and hence must be interpreted as 
a radiation reaction force. Since it generalizes Casimir's result for a situation where (at least) one of the mirrors is 
moving, it has been called motional Casimir force in Ref . |^] , where a one-dimensional calculation is presented for the 
case of partially-transmitting mirrors. We consider again the case where one of the mirrors is at rest, and then write 
the Fourier transformed force <5-F[a;] as: 

5F[<J\=x[f>]&lW (52) 

As discussed in Ref. Jg) , linear response theory provides a result for the imaginary part of the susceptibility function 
x[cl>], which corresponds to the dissipative component of thc force, in terms of thc function C ff [lü] representing the 
spectrum of fluctuations of the force operator on a standing mirror: 

ImxH = ^(CffM - C FF [-u]). (53) 

The spectrum of fluctuations CffM is defined as the Fourier transform of thc force correlation function. It may be 
written in terms of the two-photon matrix elements obtained above as follows S: 

C ff [w]=2k K"-"í-"í)\{nlj,n't'3'\F\Q)\ 2 , (54) 

{nlj.n'l'j'} 

where, as in Eq. (|3^), each pair {n£j, n'l 1 j'} is included only once (regardless of the ordering). 

The matrix elements of the force being given by Eqs. (fÏ3"l), ([Ï6| ) and (|47|), we replace the r.-h.-s. of Eq. ( |54"|) into 
Eq. (O) to find 



(f)2(^L)2 +(( A n)2+^) í 



[S(u 



<S n )-8(u, + uí+uí) . (55) 



Eq. ( |55|) p rovides the result for the dissipative component of the force exerted on the mirror. The term with n = 
in Eq. ( |Ü5| ) is particularly interesting because it allows for a comparison with the results obtained in Ref. || for a 
1D sealar field. As discussed in Appendix C, we find that in this case the two polarizations (represented by the two 
terms in the r.-h.-s. of Eq. (|55|)) give identical contributions to the dissipative susceptibility, which are in agreement 
with the perfectly-reflecting Hmit of the 1D susceptibility function derived in Ref. 

As mentioned before, I myM is directlyrelated to the number of radiated photons by energy conservation. Indeed, 
comparing Eqs. ( j3l"l ) and ( p^ ) with Eq. (p5|), we find 

£ TuoíiNZT + KJ>) = í %(lmx[u])mu}\\ (56) 



Eq. (56) shows that the energy supplied to the field by the radiation pressure force ÓF[uj], given by its r.-h.-s., is equal 
to the total radiated energy. In the next section, we discuss in detail the properties of the radiation by taking the 
specific example of sinusoidal motion. 
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IV. PHOTON PRODUCTION RATES 



In this section, we discuss in some detail the properties of thc radiation cmitted inside the cavity, starting from 
the expressions for the two-photon probabilities given by Eqs. ( pÏ8| ) and ([Ï9j), which were shown to agree with the 
results for the photon numbers N n ^ obtained directly from the moving boundary conditions and given by Eqs. ( |3Ï| ) 
and (|32|). We assume that the second cavity mirror is at rest at x — L (hence 5q2 = 0), and that the first mirror 
oscillates around x = according to the law: 

5q{t) = Sq e- W/At cos{uj t), (57) 

where the amplitude Sqo and frequency lüq satisfy the non-relativistic condition LüoSqo <C 1. Moreover, we assume that 
the damping time Aí is much larger than the period of the mechanical oscillation: 

w Aí > 1. 



We first consider thc 1D limit of the results found in sections 2 and 3, by picking up thc photon pairs with n = 0. 
Refs. fïl]l and fL2f presented a 1D nonperturbative treatment for the situation where the mechanical frequency lüq 
satisfies the resonance condition 

ujq = , (58) 

for two longitudinal cavity modes £ and £' (Ref. [Q considered the particular case £ = £' = 1, whereas Ref. |ÏÏ]] also 
considered the case £ = 2, £' = 1). We may discuss the relation between such formalism and the one presented in this 
paper by taking the Fourier transform of Eq. ( |57j ) and computing the two-photon probabilities in the resonant case 
(we omit explicit reference to polarization while discussing the 1D limit). As shown in Appendix C, we find 

TT 2 £f 

\c { oe,oi'}\ 2 = { i + 5u ,) L 2 ^ 2 ^ 2 - ( 59 ) 

According to Eqs. ( |5Ïj ) and (|59|), thc numbcr of photons grows quadratically in time in this case. The same timc 
dependence may be obtained as the short time limit of the 1D nonperturbative results found in Refs. jTÏJ and [fÏ2|| . 
Such behavior is related to the property that the spectrum of a 1D perfect cavity is discrete, and it was also obtained 
in the model of a 3D perfect closed cavity system discussed in Ref. |Ï3| . In the case of a continuous spectrum, on the 
other hand, the emission probabilities grow linearly in time as long as the perturbative approximation is vàlid, which 
is well-known from the derivation of Fermi's golden rule, so that in the end the meaningful physical quantities are 
the photon production rates, as we show below. That is the case of a partially-transmitting cavity, even in the 1D 
approximation (see Ref. 0), as well as of a 3D open cavity, as for instance the two parallel infinite plates considered 
in this paper, even under the assumption (considered in this paper) of perfect reflectiveness. 

In the 3D case, we have to sum over all possible vàlues of k|| n in order to compute the probability SP\ ^ for emission 
of a pair of photons with indices £\ and I2 and polarization j. Since ky is actually a continuous variable, we replace 



E 



s 



d 2 k\\ 



(27T) 

The probabilities do not dependent on the direction of k|| , hence we find, first for TE polarization, 

S f°° 

SP ïut 2 = 2^l dww|c{^TE,£ 2 T E }(w)| 2 , (60) 



where |c{^ lTEi £ 2TE }(a;)| 2 is obtained from Eq. J4S|): 



, _ 1 ( £ïk\ f £ 2 tt\ \Sq[u> + Co íít 



where 







represents the frequency of the 'twin' photon of index £2 and which is emitted simultaneously with the photon of 
frequency u> and index i\. We perforin the integral in Eq. (|6^) in the limit of very large Aí, so that &}[uj] is sharply 
peaked around ojq. In this case, each pair £1,^2 determines completely the freqüències lü\ and W2 of the two photons. 
Moreover, it also implies well defined vàlues for the angles between the direction of emission and the x direction, 
which we denote as 8\ and 82. In fact, we have 

ui+Lü 2 = ^Oj (62) 

as in the problem of parametric amplification by a nonlinear medium; 

u>\ sin#i = lü2 sin^2 (63) 

expressing the plane symmetry of the cavity, and which is loosely analogous to the phase matching condition in 
nonlinear òptics; and finally two additional equations which result from the boundary conditions on the two cavity 
mirrors: 

£iir 

uiiCos9i = -j-, (64) 



^1 = TT 1 + m ' ( 65 ) 



with i = 1,2. Eqs. (|62|)-(|64|) may be solved for üj\, LJ2, 8\ and 62 as functions of lüq, £1 and £2- We find 

2 l /3 2 
where 

(3 = uj q L/tt 

is the ratio between the cavity round trip time-of-flight and the mechanical period. Accordingly, the spectrum of 
photon emission, which is continuous in the case of a single moving mirror fïs|| , becomes discrete as a consequence of 
the two additional conditions, given by Eq. (|64|), and which are associated to the presence of the second mirror that 
constitute the cavity. For a given value of f3, the set of emitted freqüències is obtained from Eq. (^5|) by taking all 
positive integer vàlues of i\ and £2 in the range defined by 

h+£2<P- (66) 

In the case of TM polarization, the vàlues £\ = and £2 — are also allowed — they correspond to travelling wave 
modes propagating parallel to the plane of the mirror (waveguide modes). As for the spatial direction of emission, the 
photons are emitted along directions defined by a set of cones (whose axis of symmetry is the ï-direction) , each pair 
£1, £2 defining allowed vàlues for 8\ and 62 according to Eqs. (|62|)-(|64|). As an example, consider the value (3 = 2^/2. 
Since [3 < 3, the only allowed vàlues for TE polarization are l\ = íi = 1, corresponding to a pair with lú\ = 102 = wo/2, 
and 9\ = 82 — 45°, which is however not emitted in the case of rigid motion of the cavity. For TM polarization, on the 
other hand, there are two additional pairs: one with í\ = 1, 1% = (rigid motion), giving lj\ = 9wo/16, lü2 — 7wo/16, 
81 w 51°, and 9 2 — 90°; the other with £\ = £ 2 = (elongation motion), giving wi = ui2 — ^o/2, and 6\ = 62 = 90°. 

We compute the photon production rate for emission at a given pair of allowed freqüències assuming that the 
integrand in Eq. (^) is the product of a slowly-varying function of u> with the sharply-peaked squared Fourier 
transform of Sq(t). This amounts to replace the latter by a delta function, so that from Eq. ( p7| ) we derive: 

\5q[oj + Q lll2 ]\ 2 - 5(M 2 A£^^(c - d), (67) 
2 ujo 

and noting that u> = u>i implies il>^ 2 = W2j we find the photon production rate of TE pairs with indices £\,£2 by 
replacing Eq. (p7|) into Eq. (|6l|) and performing the integral in Eq. (BOI) : 



)TE 



S_ fhnY (hvY (Sqof 



W2k = = 375 ( ^ ) l^J — • (68) 



Note that the linear time dependence found for the probability SPJ^i 2 originates from integrating over the wholc 
width of |(5ç[íí;]| 2 , instead of taking just the peak value as in the derivation of the 1D result given by Eq. j59|). For TM 
photons, we find, starting from Eq. (E9I) and following the same method, 
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W 'ui,-ü? {1 + Seio){1 + Sl3Q) ^ ■ (69) 

We may also derive the total production rate for a given value of (3 by adding over all vàlues of í\ and £2 in the 
range defined by Eq. J66|): 

E W /U> (70) 

with j = TE, TM. In the figure, we plot W™ and W TE , both divided by the total production rate of TE photons in 
the case of a single moving mirror (see Ref. |ï^] ) , 

W™ glc = j^Si&io) 2 ^, (71) 

as functions of (3. The curves displayed in the figure for TE and TM polarizations are similar to those representing 
the decay rate of a classical dipole at the midpoint between two perfect plane mirrors along the direction parallel and 
perpendicular to the mirrors, respectively. Underlying both effects are the properties of the vacuum field in the case 
of a plane cavity geometry and the corresponding mode spectral density function |2(Í[| . The most striking differences 
between the two problems are related to the two-photon nature of the quantum radiation effect considered in this 
paper (that explains, for instance, why, as displayed in the figure, the TE photon production rate vanishes for f3 < 2, 
whereas the parallel dipole decay rate vanishes for (3 < 1 only). 

As in the problem of a decaying dipole, the photon production rates jump at integer vàlues of (3. This originates 



from adding the contribution of a new pair £1,^2 within the range defined by Eq. (66). The jumps for TE polarization 
are comparatively larger, which may be understood from the fact, discussed in detail in Ref. p8[ in the case of a 
single mirror, that TE photons are preferably emitted near the x direction, thus being more sensitive to the discrete 
nature of the wavevector along that direction. For both polarizations, the jumps become smaller as (3 increases, and 
then the curves approach their asymptotic vàlues for (3 — > 00 which are indicated by the dashed lines in the figure. As 
expected, they correspond to the photon production rates for a single moving mirror — the rate for TM polarization 
being 11 times larger than the rate for TE polarization, given by Eq. (|7l|). Alternatively, the asymptotic limits may 



be derived directly from the analytical results given by Eqs. (68) and (|69j) if we replace the sum in Eq. ( |70[ ) by an 
integral. In fact, performing the integral in the case of TE polarization leads to the expression given by Eq. (|7Ï|), 
whereas the result for TM polarization comes with an extra factor of 11. 

Of special interest is the behavior of the TM photon production rate in the range < (3 < 1, where, according to 
the figure, W™ increases strongly as [3 decreases to zero. The precise dependence on (3 may be obtained by replacing 
lüi = lü2 = üJq/2 in Eqs. (69) and ( |70| ) and comparing with Eq. (ffï|): 

45 W TE , 

W™ = W™ = -j^T 1 - ( 72 ) 
Such dependence with (3 suggests that the most favorable orders of magnitude occur for (3 < 1. In this range the 



photons have frequency lüq/2 and propagate along directions parallel to the mirrors. Following Rcfs. |2j and |18| we 
rewrite the photon production rate given by Eq. ( |72|) as: 

"""-3*( S ?)> < 73) 

where w max = cj &7o ^ s the maximum value of the velocity, and Ao = 2ttc/ujq is half the value of the wavelength 
of the emitted photons (we have reintroduced the speed of light c). As in Ref. we take v max /c = 10~ 7 and 
lüq = 2tt x 10 10 sec _1 , yielding Ao = 3cm. A real experiment would hardly employ moving mirrors with transverse 
dimensions larger than that, thus we take S/Xq « 1 in order to have a crude estimate of the orders of magnitude, even 
though diffraction effects at the borders of the mirrors, not taken into account in this paper, are of course relevant in 
this range. Finally, we take L = Ifim, giving (3 » 10~ 4 . Eq. ( f73| ) then yields W™ «4x 10 3 photons/sec. 

V. CONCLUSION 

We have calculated the photon production rates for a plane cavity with moving mirrors by two different methods. In 
the first approach, we consider the boundary conditions for perfectly reflecting moving mirrors in the long wavelength 



11 



approximation and assuming the ficld modification due to the motion to be small. We then obtain an input-output 
transformation for the field bosonic operators which allows us to compute the number of emitted photons. In the 
second approach, we start from an effective perturbative Hamiltonian and apply usual first-order perturbation theory. 
This method is considerably simpler since the expressions for the fields scattered by a moving mirror are not required, 
and establishcs a clear connection between the photon emission effect and vacuum radiation pressure. Furthermore, it 
cxplicitly unveils the fact that the photons are emitted in pairs (that satisfy simple properties expressing the symmetry 
of the plane geometry), essentially because the effect is contained in the time evolution of the field state vectors rather 
than in the evolution of the field operators. The two methods provide the same results for the photon production 
rates, hence justifying the somewhat heuristic Hamiltonian approach. 

Radiation is generated even when the distance between the mirrors is kept constant, showing that the effect is not 
simply a consequence of modulating the optical cavity length. When the initial cavity length L is much smaller than 
2irc/u)Q (we have considcrcd in detail the example of a quasi-sinusoidal motion at frequency u>o), however, radiation 
is emitted only in the case of relative motion of the mirrors, and the generation rate is enhanced as L decreases. In 
this regime, the photons are generated at the subharmonic frequency Wo/2, propagate parallel to the plane of the 
mirror, and are TM polarized. Such enhancement effect is closely rclatcd to the properties of the radiation emitted by 
a single mirror in free space fïs|] , whose spectrum for TM polarization is sharply peaked around the frequency lüq/2. 

The orders of magnitude for the photon production rate found in this paper suggest that the motion induced 
quantum radiation effect may be observed under certain conditions. However, a careful analysis of the diffraction 
cffccts ncar the border of the mirror would be necessary if a quantitative comparison with experimental results is 
required, sinec the field wavelengths involved would probably be of the order of the transverse dimensions of the 
mirrors. 
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APPENDIX A: BOUNDARY CONDITIONS FOR A PERFECT MOVING MIRROR 



In this appendix, we derive the boundary conditions in the case of a perfect plane mirror moving along its normal 
direction. We take a Lorentz frame S"(ío) whose trajectory in the laboratory frame S is given by x = Sq(to)(t — to) + 
Sq(to), so that S"(ío) represents the instantaneously co- moving frame at time to- Quantities measured in S'(to) are 
denoted by primed letters. The space-time coordinates in S'(t ) are related to those in S by 



X = 7 (x' + Sq(t )t') + Sq(t ), rj| = T|| , t = 7 (*' + Sq(t )x') + í , 

where 7 = [l — (Sq(to)) 2 ~\ • The electromagnètic fields E' and B' satisfy the following conditions: 

X x E'(.t' = 0,r||,í' = 0) = 0; x • B'(x' = 0,r\ { ,t' = 0) = 0. 
In the case of TE polarization, the condition for the elèctric field yields 

a t ,A <TE) '( a; ' = o,r; l· í' = o) = o, 



and since x ■ A (TE)/ = 0, we have from Eqs. ( |AÏ] ) and (A3) 

7 (Sq(to)d x + d t ) A (TE) (x = Sq(t ), r,, , t = t ) = jd t A (TE) (x = Sq(t = í ), r,,, t = t„) - 0. 



(Al) 



(A2) 



(A3) 



(A4) 



where dt represents the total time derivative. Since ío is arbitrary, Eq. (A4) implies that A (TE) (x = Sq(t),ru,t) must 
assume a constant value, which is taken to be zero as in Eq. (||). 



As for TM polarization, the condition on the elèctric field given by Eq. (A2) jointly with Eq. (g) yield 

d x ,A™'(x' =0,r[,t' = 0) = 0, 



(A5) 



On the other hand, we may write the l.-h.-s. of Eq. (A5) in terms of unprimed quantities by using again Eq. (Al) 
and the fact that x ■ A {TM) ' = : 



d x> A™\x' = 0, rí, , f = 0) = 7 (d x + Sq(to)d t ) A (TM} (x = &t(to),T\\ ,t = t ), 



(A6) 



and then we obtain the boundary condition as given by Eq. (|5t) from Eqs. (A5) and (A6). 



APPENDIX B: TWO MOVING MIRRORS 



In this Appendix, we consider the more general case where both mirrors move along the x direction. The first 
mirror is at x = 5qi(t), whereas the second one is at x = L + 5q2(t). As before, L represents the initial cavity length. 
For TE polarization, the boundary condition at the second mirror now reads 



A< TE '(L + <5( ?2 (í),r||,í)=0, 



(Bl) 



which yields, in the long wavelength and perturbative approximations, the following additional boundary condition 
for the motion induced perturbation i5A (TE> : 



ÍA (TE) (L,r||,í) = -Sq 2 (t)d x A^\x = L,r h t) 



(B2) 



Working in the mixed reciprocal space and using the normal mode decomposition of A^ E) as given by Eq. (|ÏJ), Eq. 
(p2h leads to 



/ P 



e=i 



ujÍSL 



8q x [u - ufàoiT 



(B3) 



Of particular interest in Eq. (B3) is the factor cos(£7r) = (— 1) that comes from evaluating the x derivative of A^ 



at x 



L. Eq. (B2) jointly with Eq. (|10| ) define a boundary value problem to be solved with the aid of the Green 



13 



functions given by Eq. (|23). We first employ the retarded Green function G^'J i (x, x') to solve for the total field AJ™ 1 
in terms of the input field AÍ™^ : 

A< TE, M = A^Ocw) +6A^ ) (L,u > )d x ,G^ R (x' = L,x;u>) 

-8A^\Q^)d x .G^ R {x' = 0,x;u). (B4) 



As in Sec. 2, we also solve Eq. (|B3j) in terms of the output field A { out ' n with the aid of the advanced Green function 
G^iy(x, x 1 ). The connection between output and input fields is then provided by the difference 

dx'Gn^(x' = L,x) — dx'G^lj (x' — L,x) 

= -?B- 1 ) í (|) sin (T)?( í(w - w " ) - í(w+ ^· (B5) 



As explained in Sec. 2, we derive the linear transformation between output and input TE bosonic operators from 

E q s. PD-dr 



oo 

(TE) (TE) * 



E 



^tt\ /f7r\ 1 



{(<%[^-c4']-(-l)^ : 



'2 K,, 



'in, 



(íç^ +<4\ - {-If+e&totà + utí) « E l^) + } . (B6) 

From Eq. ( fB6| ) we may calculate the number of photons iV^™' by taking the average of the output number operator 
over the input vacuum state as in Eq. (^0|). For TM polarization, we extend the method employed in Sec. 2 to take 
into account the motion of the second mirror exactly as discussed above for TE polarization. 

Alternatively, we may compute the photon numbers from the effective perturbation Hamiltonian 

6H = -5q 1 {t)F 1 -5q 2 {t)F 2 , (B7) 

where Fi is the force exerted on mirror i by the vacuum field. Following the procedure outlined in Sec. 3, we derive the 
two-photon creation probabilities given by Eqs. ( [48| ) and (^9|). As in the case of a single moving mirror, the results 
obtained through this method arc in full agrccmcnt with those obtaincd directly from the boundary conditions. 

APPENDIX C: PHOTONS EMITTED ALONG THE NORMAL DIRECTION 

In this appendix, we consider in detail the contribution of the degenerate two-photon states in the derivation of 
the photon numbers and of the susceptibility function. First note that degenerate two-photon states necessarily 
correspond to propagation along the direction perpendicular to the plane of the mirror, i.e., they are of the form 
\n = l j, n = £ j). The degenerate two-photon matrix elements of the force o per ator are calculated from the 
representation of the force operator in terms of the vector potentials, given by Eq. (|42j), and from the normal mode 
decompositions given by Eqs. (Tl3) and (|Ï3|): 



lirh 

(n = £ te, n — ^te|F|0) = -{n = £ tm, n = £tm\F\0) = -= — . (Cl) 

y/2L 2 

These results are smaller than the vàlues of the expressions given by Eqs. ( fÏ6| ) and (^7|) at n = and £ — £' by a 
factor of y/2. From them, we easily compute the degenerate two-photon probabilities by using Eq. (|40|), allowing us 
to write the correct expression for n = : 

12 ,2 1 i c A e7r ï 12 frUí y 

|C{ £te,0€'te}I - |C{oítm,0£'tm}I - 7^ ~ §^A£2 l T/ J \~T~ J I ^T/ + ~£~J I " ^ > 
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According to Eq. (|C2|), the degenerate two-photon probabilities are one-half the value found when replacing the 
vàlues n = 0, l = i' , and Sq 2 = in Eqs. © and ©. 

The contribution of degenerate two-photon states is found from Eq. (pQ): 



N J - 



E \ c {otjM'j}\ 



2\ c {0tj,0tj}\ 



(C3) 



The factor two multiplying the degenerate two-photon amplitude in the r.-h.-s. of Eq. (C3) cancels the additional 
factor one-half appearing in Eq. (02) for £ — £' , then yielding a result in full agreement with Eqs. ( [jï|) and §ï%). 
For the specific example of motion given by Eq. (p7[), and as sum ing that the mechanical frequency uq satisfics the 
resonant condition as given by Eq. (pq), we derive from Eq. (C2) the expression for the production rate of photons 
with n — given by Eq. ( |59| ) 

Since the results for the photon numbers are not modified when taking into account the degenerate two-photon 
states, we expect that the formula for the dissipative component of the susceptibility function, given by Eq. (|5E|), 
should also be vàlid for n — 0, so as to preserve the connection between dissipation and total radiated displayed by 
Eq. ( |56| ). In fact, we may write separately the contribution of degenerate two-photon states to the sum over pairs 
{nl,n'£'} in Eq. (||): 



Cff M = ^E E S(iü-Lü e n -^)\{níj,ne'j\F\0 
j n,e,e' 

+ 2?r E E % - 2£ir/L)\(0£j, 0£j\F\0}\ 2 , 



(C4) 



3 l 



where ^* t e , represents the sum over all possible vàlues of n, £ and £' excluding those where simultaneously n = 
and £ — £' . As before, the factor one-half found for the degene rate two-photon matrix element (given by Eq. (|C2^)) 
is cancelled by the factor two appearing in the r.-h.-s. of Eq. (C4). Hence we may write the expression for the 1D 
dissipative susceptibility function, Imxiofa;], by selecting directly from Eq. (p5[) the terms with n — : 



7T 3 h 



ImXioM = ^ E E U> ^ ~ + £ >/ L ï - S ^ + ( £ + e >l L ^ 



(C5) 



i=i í'=i 



As for Ref . P] , the result for the (complete) susceptibility function in the perfectly reflecting limit and in the particular 
case where only one mirror moves reads 



XM = t~ 



ILÚ 



(z) Hirr^ 



(C6) 



In order to compare Eqs. ( |C5| ) and (| 

Imx[w] = 



we must take the imaginary part of xMj then yielding, after some àlgebra, 



h 








67T 


. 2 


h 5 


(i) 



?i(rt 2 — l)d(u> — niï jV) 



(C7) 



The first term in the r.-h.-s. of Eq. (C7) represents the contribution of the ficld outside the cavity, being equal to half 
the value found for a two-sided single mirror in onc-dimensional vacuum Q. The seco nd ter m, on the other hand, 
represents the contribution of the intracavity field, which is, by inspection of Eqs. ( |C5| ) and (C7), equal to half the 
value found from taking the 1D limit in Eq. (|55|), the factor of two being related to the two polarizations taken into 
account in the electromagnètic case. 
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Figure Caption 



Total production rates of TE (a) and TM (b) photons as functions of (5 = lüoL/tt, which represents the ratio between 
the round-trip time of flight and the mechanical period. The scale of the vertical axis is such that the value one 
corresponds to the generation rate of TE photons for a single moving mirror. The dashed lines provide the asymptotic 
limits for large (3. They show that the photon production rates approach the vàlues corresponding to the case of a 
single moving mirror in this limit (note that the single mirror TM photon generation rate is larger than the single 
mirror TE rate by a factor of 11). 
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